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Abstract 

We study the non restoration of symmetries with a local order parameter in field theory at 
finite temperature. After giving an interpretation of the phenomenon, we show that hierarchy 
problems are a necessary condition for its realization in renormalizable theories. We then use 
a large N treatment, and find that high temperature symmetry can stay broken in this limit 
(in opposition with a previous result), and further that the running of couplings reinforces the 
effect in the simplest model with two scalars. 
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1 Introduction 



Since soft pion theorems in the sixties and renormalization of the standard model in the sev- 
enties, spontaneous symmetry breaking has become an inevitable building block in particle 
physics models to understand the small mass of certain scalar particles, or the large mass of 
certain gauge ones. It was later noted [|], [| that in analogy with ferro-magnetism, heating up 
such models tends to restore the broken symmetry. There are two related intuitive arguments 
why this should be expected. The first one appeals to thermodynamics: for higher tempera- 
tures, the minimum of the free energy becomes determined less by energy (which has broken 
symmetry minima), than by entropy (which intuitively is maximum in symmetric configura- 
tions). Alternatively, one can say that thermal fluctuations are able to cross the potential 
barriers surrounding the broken minima, and feel the global symmetry of the theory. 

This implies that any symmetry broken today was once unbroken when the cosmological 
temperature was high enough, and there must have been a phase transition where this freezing 
of symmetry occurs. During this freezing, various topological defects can be created. Certain 
of these, like monopoles or domain walls are experimentally excluded, which gives constraints 
on the particle physics models. For instance, this makes it hard to explain the breaking of 
discrete symmetries (like CP) dynamically: at the phase transition where the symmetry gets 
spontaneously broken, domain walls separating the regions with different discrete minima would 
be formed, and these could quickly dominate the energy density of the universe. 

To relax these constraints, it has been noted || that domains would not form if the symmetry 
remains broken, however high the temperature. This is symmetry non-restoration (SNR). 

In one of the founding papers of finite temperature field theory, Weinberg already noted 
that high temperature symmetry restoration was not an unavoidable fate in field theory, and 
produced the following simplest example. Consider a scalar theory with 2 fields in the vector 
representations of an O(Ni) ®0(N 2 ) global symmetry group. The most general renormalizable 
potential is then: 



Notice this does not exclude a negative mixed coupling A12, provided it is not excessively large. 

For high enough temperatures and small enough couplings, the effective mass corrections 
are dominated by the quadratic divergence of the one loop tadpoles: 

m 2 2 (T) = ™2 + ^ [(N 2 + 2)A 2 + iViAia] (1 + 0(A, m/T)) (3) 
and similarly for (1 <-> 2). Thus, for 



V{4>) = \m l 2 <p 1 2 + \m 2 2 <p 2 2 + ^Ai0i 4 + -\ 2 <p 2 4 + lXi 2 <Pi 2 <p 2 2 , 



(1) 



for which boundedness at the classical level requires 

A1A2 > \i 2 , Ai, A2 > 0. 



(2) 




(4) 
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any increase in temperature tends to increase the thermal expectation value of 02 



<*> = ~ T (5) 

up to arbitrarily large values, while the potential stays bounded from below. 

To get an intuitive understanding of this disturbing result, it helps to plot the equipotentials 
in the (0i,02) plane. To simplify, we will consider the limit of large fields relevant for high 
temperatures, and consequently neglect the mass terms in this plot. The thin dashed line is an 




Figure 1: Contour plots of the classical potentials (0), for A12 = (thin dotted line) and A12 
close to — \i2™ x (thin plain line). Thick grey disks represent thermal fluctuations about the 
cold (dotted line) or hot (plain line) vacuum expectation value. 



equipotential for the case of vanishing mixed coupling A12. Turning on a negative A12, we see 
the thin plain equipotential grows diagonal extensions that would eventually make the potential 
bottomless for — A12 > A^ a2 . Since all other equipotentials are homothetical to the ones drawn, 
the minimum of this classical potential lies at the origin. However at finite T, we should not 
only minimize the potential energy for a point, but also for the thermal fluctuations around 
it, roughly depicted as thick grey disks of size ~ T. Because of the pinch in the vertical 0i 
direction, the dashed oval get pushed away from the origin, to reach the finite T expectation 
value surrounded by the plain oval. 

With this picture in mind, symmetry non restoration seems to be an indubitable property 
of certain scalar theories, at least when couplings are small enough that one loop calculations 
can be trusted. It is thus not surprising that certain crystals (which typically have a rich 
scalar, vector, tensor, ... excitation spectrum) undergo such an inverse transition where the 
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higher temperature lattice has less symmetry than the lower T one. However, this comparison 
also calls for a distinction between inverse symmetry breaking and symmetry non restoration. 
Indeed, by increasing temperature enough, these crystals inevitably melt, thus restoring the full 
spatial group in the liquid phase. This stresses that any claim for symmetry non restoration is 
a strong statement about the ultra-violet structure of the theory. 

In section |2], we will see that for inverse symmetry breaking to be possible in a renormal- 
izable theory, its scalar sector must necessarily suffer from hierarchy problems. This excludes 
renormalizable supersymmetric theories. 

In section |3|, we try to see whether there is any sign for symmetry restoration when the 
model (|1|) is heated close to the point where its couplings get uncontrollably large (as must 
somehow happen for crystal excitation modes close to melting). We find contrarywise that the 
running of couplings seems to increase the symmetry non restoration region. 

2 The U.V. Price for Inverse Symmetry Breaking: Hi- 
erarchy Problems 

Consider the most general renormalizable field theory, containing scalars (pi, fermions ip a and 
vectors A a . It is easy to prove that if the theory is hierarchy safe, it cannot undergo inverse 
transitions. By "hierarchy safe" we here mean the naivest possible definition, namely that 
relations between dimensionless couplings prevent all dimensional parameters in the lagrangian 
from receiving power-like divergencies at one loop. 

Let us in a first step exclude scalar singlets. To decide whether symmetry is broken or not, 
we must study the stability of = which is determined by the lowest non-trivial terms of 
the potential. Although scalar trilinear couplings exist which could give rise to linear terms at 
one loop, group theory forces such linear terms to vanish. We must thus focus on mass terms, 
which by the above hierarchy requirement only receive logarithmic radiative corrections. The 
quadratic divergencies exactly cancel at one loop between bosons and fermions: 

A 2 

= (c ijb - c ij)^^ + O(logAtfy); c ijb = c iJ} (6) 

where c^&m is a combination of all couplings of scalar i to bosons (fermions). Such a cancellation 
is for instance automatic in supersymmetric theories f|, but the argument here is simpler and 
more general. 

Turning on temperature, we find 

m 2 (T)-m 2 (0) = c ^ + ClJ Z! + 0(logT) (7) 

= 3cv^ + 0(logT), (8) 

using "hierarchy safety" to obtain the last line. Contrary to bosonic contributions like ([3]), Qj 
can easily be checked to be positive semi-definite. Indeed, they are the eigenvalues of a mass 
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matrix of the form: 

9ij = yiabVjab = J2 ViAVjA' (9) 
a,b A=(a,b) 

which can be seen as the metric induced on the space of bosonic fields i, from the identity 
metric 5ab in the space of fermionic pairs A, by the linear embedding provided by Yukawa 
couplings DiA- By suitable rotations in the space of bosonic fields, this hermitian matrix can 
always be diagonalized: 

9ij = $ij9i = $ij \ViA\ 2 ~ dijCij- (10) 

A 

This proves that c,ij are positive as claimed, and furthermore, that they can only vanish for 
scalars totally decoupled from all fermions. Thus, ((pi) = is a local minimum, whose stability 
increases as temperature is raised. 

Let us now consider scalar singlets, which allow for new types of linear and trilinear vertices. 
The only way these can alter the previous conclusion is by having some singlets developping 
a quadratically divergent expectation value which could contribute to (|B|) through trilinear 
coupling. However, this is forbidden by our naive definition of hierarchy safety as it induces 
an even stronger problem than the usual one for scalar masses. One can for instance show 
that certain fermions then necessarily receive quadratically divergent contributions to their 
mass. Hierarchy problems in the scalar sector are thus necessary conditions to have symmetry 
no n- restoration, or even just inverse symmetry breaking. 

3 Gap Equations and Running Couplings 

Let us now focus back on Weinberg's simple model ([!]) and consider the question of coupling 
running. Neglecting mass terms, the zero temperature one-loop renormalization of the couplings 
reads (dt = d\og k uv / IQn 2 ) 

(N 1 + 8)X 1 (t) 2 + N 2 X 12 (t) 2 

(N 2 + 8)X 2 (t) 2 + N 1 X 12 (t) 2 (11) 

X 12 (t)[(N 1 + 2)Ai(*) + (N 2 + 2)A 2 (t) + 4A 12 (t)] 

Independently of the sign of Ai 2 , all couplings thus tend to grow stronger with increasing 
Auv, just like for a single self-coupled scalar. As implied by triviality in the latter case, the 
theory is only defined below A max , some physical maximum UV cut-off, where the couplings 
start blowing up. 

The presence of such a cut-off requires some care in the definition of symmetry non- 
restoration. Indeed, the existence of a maximal energy scale means that symmetry (non-) 



d\i 
~dT 
d\2 
dt 
dXu 

~df 
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restoration can only be probed up to that temperature, and no easy claim can be done about 
what happens beyond, where the theory becomes intractable?]. Furthermore, an upper cut-off 
opens new possibilities: taking for instance a single scalar with — m 2 (T = 0) ^> XK max 2 will 
forbid symmetry restoration up to T < A max . This trivial effect is however qualitatively differ- 
ent from what we have with (H), where a scalar thermal expectation value constantly increases 
with temperature. 

We will thus focus our attention on the limit T ^> m(T = 0), for which we may as well 
take vanishing masses to start with, and will call symmetry non restoration in the presence of 
a UV cut-off, the possibility that the symmetry breaking thermal expectation value shows no 
tendency to decrease as temperature is increased all the way up to the cut-off. In the simplest 
approach @, this happens provided a condition on the couplings is satisfied: 

C{Xi, A 2 , Aw) = (iV 2 + 2)A 2 + jViA u < 0. (12) 

However, if couplings are larger, the thermal masses (Q) get comparable with T. To keep 
properly track of 0(m/T) terms, (|3]) can be replaced, in the large (Ni,N 2 ) limit, by self- 
consistent gap equations 0, also known as bubble, super-daisy, or cactus resummation: 

^vm = i^ Xlf(Xl)+ ^ f(X2) 

A ±«m = ^±v X2f{X2)+ ^ fiXi)i (13) 

/<*) 4/" -^(e^ - I)- + OC) (14) 

It JO VP + x 7r 

is exponentially decreasing for large x. The approximation / = 1 just gives back (§), while 
expanding f(x) around the origin gives perturbative corrections for any N 0, [5|. The new 
region of symmetry non restoration is now bounded by x 2 = which, upon using (|13 ) and 
dropping 1/N corrections, becomes: 

C(X h A 2 , A 12 ) =\ 1 -2*+ —^f-^-N^/N^) 2 < 0. (15) 

A 2 iV 2 A 2 

This condition on the couplings contrasts with [0] , which claims the SNR region to be empty 
in the same large iV limit. The origin of this difference is that we kept all terms in ([13]) that 
stay finite in the limit AT 1)2 — > oo with Ai i2j i 2 iVi )2 finite. One can further explicitly check that 
there is no x 2 > solution inside the region (|l5|) . 

Notice that, just like (0), (0) purely gives a condition on the couplings which are kept 
constant in the gap equations approach. Gap equations can not produce a critical temperature 
because they have no energy scale other than T in the large temperature limit. 

3 Notice this equally casts doubts on the restoration of the symmetry above that temperature 



where 
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Figure 2: The Symmetry Non Restoration (SNR) domain in couplings space: the light grey 
area is given by naive one-loop calculation; the dark grey region is what remains according to 
gap equations; all stable models should lie below the thick curve. 



Taking couplings evolution into account will bring in the scale A max alluded above. To 
get a feeling for the direction in which this might go, let us notice that thermal masses come 
predominantly from "hard" momenta of order T in expressions like (|14]) . It is therefore natural 
to take the couplings showing up in flT3|) as the running couplings at scale T. For a given 
theory fixed e.g. by the values of "bare" couplings at a scale as close as reasonable to A max , 
the running from that initial scale down to T is close to the running we would have at zero 
temperature. Indeed, the largest effects of temperature show up in the infra-red, for scales lower 
than T where the running effectively becomes 3 dimensional |10[ . We can therefore tentatively 
use the same equations flTTI ) and now interpret dt as d\ogT/l6ir 2 . This procedure was recently 



legitimated, at least for small couplings [11 



It is interesting to notice that the boundary of the stability region A1A2 = \\ 2 ls invariant 
under the large N running of couplings. This is reassuring about the use of the same condition 
at all scales. Indeed, the global stability of a theory cannot change by integrating out degrees 
of freedom. Hence, the finite N result that the flows of A can cross the "stability" boundary 
must mean that operators higher than quartic become relevant to this question [|EJ]. 

At every temperature, the symmetry non restoration region is thus given in terms of a 
condition on the running couplings at that temperature. We will now show that increasing T 
brings new points in this region, so that symmetry non restoration is effectively enhanced by 
couplings running. This region is bounded by the invariant stability condition and by 



Ai = A lc (A 12 , A 2 ) = ^ - -^-^(-JVaAa/iVxAu,) (16) 
where g(x) = / _1 (x) 2 . Let us compare the evolution of Ai as given by ([TT|) , with the evolution 



6 



along the surface given by: 



dXic ^ d\ lc d\ l2 OX lc dX 2 
dt ~ 8X 12 dt d\ 2 dt [ ' 

= N^l + N^ + ^-M-^-) (18) 

dX\ X 12 f N 2 X 2 

+ Tr ^gg(- Tr ^) (19) 



dt N 2 X 2 y iViA 
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The last term of this expression is positive definite, since A12 is negative, and g is purely 
decreasing just as / and As a result, when increasing t, a point on the surface Ai c (Ai2, A2) 
is moved to a point with a smaller Ai than its projection on the surface. This means that the 
domain of symmetry non restoration increases with temperature. 

Let us now critically examine the approximations used to obtain this result, in comparison 



with the exact renormalization group equations used in [11|. We have first neglected the running 



of masses as a function of the momentum scale in the gap equations. Taking this into account 
would change the function f(x) in (0), but it could not change its purely decreasing behavior, 
which was the only information needed to reach our conclusion. 

We have further neglected masses in the running of couplings (n|). These would introduce 
different factors for the loops with different masses. However, the remarkable cancellations 
leading to ( |19D would still work, leaving just an mi-depending factor in the last term of that 
equation. Since the bare theory needs a large negative m\ to achieve a small infrared mf, this 
factor is dominated in the large N limit by the N± — 1 massless Goldstone modes, for which 
(|19l) is thus unaltered. 

Finally we have kept couplings fixed in the gap equations. This is a good approximation if 
couplings evolve slower than masses, like at zero temperature where their logarithmic running 
is negligible compared with the power-like mass evolution. However this assumption clearly 
breaks down when nearing the triviality pole. It should be noted that for T ~ A max , we loose 
control over the theory anyway, as unavoidable regulator dependencies start creeping in. This 
shows up in the exact renormalization group approach [|TT|] under the form of strong oscillations 



in the T dependence when T/Auv is not small. These prevent any firm conclusion in this 
region, and the approximation made at least has the merit of cutting them. Going beyond 
our approximation requires either knowledge of the theory beyond A max , or at least a way to 
reach higher temperatures. This is possible by introducing a different spatial and temporal 
cut-off. The spatial one is still limited by the same A max , but the temporal cut-off can now 
be separately sent to infinity, thus allowing for arbitrarily large temperatures. Physically, this 
corresponds to measuring temperature by the tiny energy dependence in the distribution of 



modes much smaller than T. It has been proven on the lattice |13| that this leads to symmetry 
restoration. Although this is far from the continuum limit, we expect a similar conclusion with 
a momentum cut-off. 

Outside of this admittedly artificial regularization procedure, the theory can only be defined 
beyond A max by introducing new physics which may restore the symmetry, or not. The only 
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known fact is that in the latter case, the hierarchy problems get worse beyond A max . 

4 Conclusion 

In this paper we studied the non restoration of symmetries (SNR) with local order parameters. 
We showed that in the framework of renormalizable theories, inverse symmetry breaking (and 
thus a fortiori SNR) necessarily rested on a hierarchy problem for the scalar field associated 
with the broken order parameter. This excludes renormalizable supersymmetric theories. It 



has recently been shown |TJ[ that including non-renormalizable operators does not remove this 
constraint. 

We have then shown that in the simplest model of Weinberg ([I]), the inclusion of running 
couplings does not lead to symmetry restoration, even when venturing into the ultraviolet strong 
coupling regime by means of large N techniques. On the contrary, the running of couplings 
tends to enlarge the region of SNR, in accordance with the recent result of |J on the effect of 
next to leading order corrections, the first ones sensitive to coupling running. 

This all should not obscure the fact that this model possesses an intrinsic maximum energy 
scale, which translates into a maximum temperature. Such a maximum scale is quite general 
in theories with scalars, and is for instance unavoidable when the scalars fully break all non- 



Abelian gauge symmetries |]T5j . Without pretending to too much generality, we can say we failed 
to construct a model with SNR where all couplings enjoy asymptotic freedom. The difficulty 
lies in that the scalar couplings must be large enough to dominate certain scalar self-energies. 

This maximum scale poses a problem for the application of SNR to eliminate topological 
defects in cosmology. Since classical cosmology should be defined for temperatures up to the 
Planck scale, we have the following dilemma. Either we tune parameters so that the maximum 
scale is pushed above the Planck one. In that case, it could be argued that the classical universe 
is born in the broken phase, but one should somehow cope with graviton loops, and show they 
do not alter this statement. Or SNR only holds up to some temperature below the Planck 
scale. In that case, inflation is needed to dispose of the topological defects created at that 
temperature, and the role of SNR is just to push the moment where these defects appear back 
in time, earlier than inflation. However, the higher in energy we push them up, the more 
critical the necessary hierarchy problem becomes. This seems an unavoidable ultra-violet price 
for SNR. 

To conclude, let us stress that we did not touch here symmetry realizations with a non- 
local order parameter. For these, the picture is totally different. For instance, a duality 
transformation can exchange high and low temperatures. However, the order parameter broken 
in the high temperature phase is then non-local. Similarly, the Polyakov loop is a non-local order 
parameter for the deconfining phase transition, and its high temperature thermal expectation 
value breaks a center symmetry, albeit in a controversial way [IB, |TT| . 
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